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Abstract 

The impact of higher-order final-state photonic corrections on the precise determination of the 
VF-boson mass at the Tevatron and LHC colliders is evaluated. In the presence of realistic selection 
criteria, the shift in the W mass from a fit to the transverse mass distribution is found to be about 
10 MeV in the W — > \xv channel and almost negligible in the W — > ev channel. The calculation, 
which is implemented in a Monte Carlo event generator for data analysis, can contribute to reduce 
the uncertainty associated to the W mass measurement at future hadron collider experiments. 
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Precision tests of the Standard Model require a more and more accurate knowledge of the 
basic parameters of the theory. In particular, future measurements of the W^-boson and top 
quark masses at the Tevatron and the LHC colliders are expected to considerably improve 
the present indirect bound on the Higgs-boson mass from electroweak precision data. As 
recently discussed in the literature a precision of 27 MeV (16 MeV) for the W mass M w 
is the target value for Run Ha (Run lib) of the Tevatron. An accuracy of 15 MeV is the 
final goal of LHC Q. 

In order to measure My/ with such a high precision in a hadron collider environment, it 
is mandatory to keep under control higher-order QCD and electroweak radiative corrections 
to the W and Z production processes. The status of QCD corrections to weak boson pro- 
duction in hadronic collisions is reviewed in Ref . jjj , while recent progress in the calculation 
of electroweak corrections, as achieved by means of independent calculations DDQ , is 
summarized in Ref. [7( . As shown in Refs. , electroweak corrections are dominated 

by photon radiation effects and, in particular, by final-state photon emission, which gives 
rise to collinear logarithms of the form a / 7tlog(s / mf) , where s is the effective centre of 
mass (cm.) energy and mi is the mass of the final-state lepton. This poses the question 
of the impact of higher-order (i.e. beyond order a) leading logarithmic corrections due to 
multi-photon radiation. A first attempt toward the inclusion of O (a 2 ) QED corrections 
was the calculation of the double-bremsstrahlung matrix elements qq' — > W — > Ivyy and 
qq — > 7, Z — > Z + Z - 77 (I = e,ji) performed in Ref. The aim of the present work is to 
evaluate the impact of higher-order final-state QED corrections on the W mass determina- 
tion at hadron colliders, by including both real bremsstrahlung and virtual corrections. To 
this end, a Parton Shower (PS) approach in QED |9j is employed to simulate multi-photon 
radiation effects. An independent calculation of multi-photon radiative corrections in lep- 
tonic W decays has appeared very recently jlo| . but without quantifying their impact on the 
W mass measurement. The uncertainty in the W mass due to higher-order QED effects is 
presently estimated by the CDF collaboration at the Tevatron to be 20 MeV in the W — > ev 
channel, and 10 MeV in the W — > \w channel 1 
by the D0 collaboration to the W — > ev channel 

From our analysis, we find that the shift in the fitted W mass is about 10 MeV in 
the W —>■ \w channel and almost negligible in the W —>■ ev channel. The inclusion of the 
present calculation in the experimental analysis would reduce significantly such a theoretical 



. An uncertainty of 12 MeV is assigned 
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uncertainty in future improved measurements of the W mass at hadron colliders. 

An appropriate theoretical tool to computephotonic radiative corrections in the leading 
log approximation is the QED PS approach |9j. It consists in a numerical solution of the 
QED Gribov-Lipatov-Altarelli-Parisi evolution equation for the charged lepton Structure 
Function D(x,Q 2 ) in the non-singlet channel. The solution can be cast in the form [9| 

D(x,Q 2 ) = U{Q 2 ,m 2 )5(l-x) 

+ j%(s,s')^U(s',m 2 ) ^ dyP(y)5(x - y) 
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where Il(si,s 2 ) = exp 



is the Sudakov form factor, P(z) is the 



e —>■ e + 7 splitting function and x + is an infrared regulator, separating the soft+virtual 
region from hard bremsstrahlung. Equation allows to compute D(x,Q 2 ) by means of 
a Monte Carlo algorithm which, as shown in Ref. simulates the emission of a shower 
of (real and virtual) photons by a charged fermion and accounts for exponentiation of soft 
photons and re-summation of collinear logarithms due to multiple hard bremsstrahlung. A 
remarkable advantage of the PS algorithm with respect to a strictly collinear approximation 
is the possibility of generating transverse momentum of fermions and photons at each 
branching, thus allowing an exclusive event generation suitable to implement experimental 
cuts according to a realistic event selection. The generation of transverse degrees of freedom 
can be performed according to different recipes, as described in detail in Ref. 0. Here, we 
generate photon angular variables according to the leading pole behavior 1/(1-/5; cost?/ 7 ), 
where fli is the lepton velocity and ^ is the relative lepton-photon angle. 

A simple recipe to evaluate final-state corrections to pp — > W ^ vl consists in attaching 
a single structure function D(x, Q 2 ) to the lepton coming from the W decay. Needless to say, 
this amounts to neglect photonic corrections due to initial-state radiation, initial-final-state 
interference and H^-boson emission. However, it is known that radiation from an internal 
off-shell particle can not contribute to leading logarithmic corrections, which are the main 
concern of the present study. Further, initial-state photon radiation requires an appropriate 
treatment, since radiation off quarks gives rise to quark mass singularities which, as discussed 
in Refs. □ DDQ , can be reabsorbed in Parton Distribution Functions (PDF), in analogy 



TABLE I: Comparison between the present calculation (HORACE) and WGRAD |g, for the 
pp ^ W —> lu, I = e, n cross sections (in pb), at the Tevatron Run II (yfs = 2 TeV) and the LHC 
(s/s = 14 TeV). 
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HORACE exponentiated 
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to gluon emission in QCD. These considerations imply that the treatment of final-state 
photon radiation alone is not gauge invariant. Nevertheless, it can be easily checked, by 
comparing the PS spectrum with the gauge-invariant factor for collinear photon emission by 
a fermion jl^ . that gauge violations are confined to the next-to-leading order corrections, 
which are beyond the approximation of the present analysis. This issue is further discussed in 
the following, by means of a quantitative comparison with the independent, gauge-invariant 
calculation of Ref . . 

In order to quantify the effect of higher-order final-state QED corrections on the W 
mass determination, we developed a Monte Carlo event generator following the approach 
described above and performed a number of Monte Carlo experiments. The technical details 
of the Monte Carlo code HORACE (Higher Order RAdiative CorrEctions) will be presented 
elsewhere. Before the phenomenological analysis, we performed a tuned comparison between 
the predictions of HORACE and those of WGRAD , to verify the accuracy of our 

calculation. The results of such a comparison are shown in Tab. |1J using default PDFs, 
input parameters and cuts as in Ref. ^J. WGRAD includes the full set of 0{ot) electroweak 
radiative corrections to W production (second line in Tab. HJ but it also gives the possibility 
to select the effect of a gauge-invariant subset due to final-state corrections (WGRAD final- 
state in Tab. nj). Therefore, the difference, at a few per mille level, between WGRAD and 
WGRAD final-state points out, when comparing with the Born predictions, the dominance 
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of final-state radiation within the full set of 0{a) electroweak corrections. On the other 
hand, it can be seen that the predictions by WGRAD final-state are in very good agreement 
with our results by HORACE 0(a), which is an order a expansion of the complete PS 
algorithm. Since the differences are well below the 0.1% level, this comparison demonstrates 
that the gauge-invariance violations present in our approach are numerically negligible. The 
contribution of higher-order effects can be seen by comparing HORACE 0(a) with our 
complete predictions given by HORACE exponentiated. Their effect on the integrated cross 
section is tiny, within 0.1%. For the sake of completeness, we performed also comparisons 
between WGRAD and HORACE at the level of differential distributions, such as lepton 
transverse momentum and transverse mass distributions, finding perfect agreement. 

Having established the physical and technical accuracy of our calculation, we move to the 
analysis of the W mass shift due to higher-order corrections. The input parameters used in 
the simulations are: 

m Ul =0 m e = 0.511 x 10" 3 GeV 

m M = 0.10565836 GeV 

a' 1 = 137.03599976 G M = 1.16639 x 10" 5 GeV" 2 
a s = 0.1185 

M w = 80.423 GeV M z = 91.1882 GeV 
•2. M 2 W _ 3G,M^ f 2a s \ 

We adopt the scheme and fixed-width scheme in our calculation. At the parton level, we 
consider the processes 

u + d -> W + -> l + + vi u + s -> W + -> l + + v x 
c + d -> W + -> l + + n c + s -> W + -> l + + ui (3) 

and their charge conjugate, with I — e,fi and CKM matrix elements according to Ref. j"!^ . 
The results for the processes pp — > W — ► I + v (Tevatron) and pp — > W — > I + v (LHC) 
are obtained by convoluting the parton-level matrix element with CTEQ6 PDFs Q]. The 
virtuality scale Q 2 is set to be Q 2 = s, s being the effective cm. energy after gluon radiation, 
in both PDFs and lepton Structure Function. The cm. energies considered are yfs = 2 TeV 
for the Tevatron and y/s = 14 TeV for the LHC. 
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To model the acceptance cuts used by the CDF and D0 collaborations in their W mass 
analyses, we impose the following transverse momentum (pt) and pseudo-rapidity (77) cuts: 

p T (/)>25GeV |77(Z)|<1.2 $t > 25 GeV (4) 

However, in order to perform a realistic phenomenological analysis and study the dependence 
of the W mass shift from detector effects, we implement, in addition to the above cuts, 
the lepton identification requirements quoted in Table I of Ref. 0]. According to these 
criteria, electron and photon four-momenta are recombined for small opening angles between 
the two particles, consistently with a calorimetric particle identification, while muons are 
identified as hits in the muon chambers with an associated track consistent with a minimum 
ionizing particle. Furthermore, we simulate uncertainties in the energy and momentum 
measurements of the charged leptons in the detector by means of a Gaussian smearing of 
the particle four-momenta, using as standard deviation values the specifications relative to 



electrons and muons for the Run II D0 detector 

The strategy followed by the CDF and D0 collaborations to extract M\y from the data 
is to perform a maximum likelihood fit to the transverse mass distribution of the final-state 
lepton pair or to the transverse momentum of the charged lepton. Here we consider the 
transverse mass, which is the preferred quantity to determine the W mass and is defined as 

M T = y/2p T {l)p T {u){l-cos<f) 1 -) (5) 

where pril) and Pt{v) are the transverse momentum of the lepton and neutrino, and <p lv is 
the angle between the lepton and the neutrino in the transverse plane. The transverse mass 
distribution, as obtained by our simulation, is shown in Fig. at y/s = 2 TeV. The distri- 
bution without lepton identification requirements and smearing effects (solid histogram) is 
compared to the distribution including lepton identification criteria (markers) and detector 
resolutions (shaded histogram). The shape of the M T spectrum is considerably modified by 
detector resolution effects, in agreement with the results shown in Refs. Q,^^. The arrows 
in Fig. n select the range 65 GeV < < 100 GeV, which is used by CDF collaboration in 
its W mass analysis and we also adopt in the fitting procedure described below. 

To evaluate the shift induced by higher-order corrections on the W mass, we perform 
binned x 2 fits and binned maximum likelihood fits to the Mj- distribution, in complete 
analogy with the experimental fitting procedure. Here we show only the results of the x 2 
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FIG. 1: The transverse mass distribution without lepton identification criteria and detector reso- 
lutions (solid histogram) , with lepton identification criteria (markers) and with detector resolutions 
(shaded histogram), in the W — > ev channel at yfs = 2 TeV. Arrows indicate the considered fit 
region. 

fits, because the results of the maximum likelihood fits are in perfect agreement with the 
former. Using HORACE, we generate a sample of pseudo-data and calculate with high 
numerical precision the my spectrum (binned into 100 bins) at the Born level and for a 
fixed, "physical" value of the W mass, i.e. = 80.423 GeV. Next, we compute the ttit 

spectrum including 0(a) leading-log corrections for 20 hypothesized W mass values, with a 
spacing of 5 MeV for the W — > ev channel and 10 MeV for the W — > fiv channel. We then 
normalize the spectra within the fit interval and we calculate, for each My/ value, the x 2 as 

X 2 = Vi,Born) 2 /{Avt a + &V 2 Born) (6) 

i 

where cr^Born an d <Ji, a are the Monte Carlo predictions for the i th bin at the Born and 0(a) 
level, respectively, and Acr^B^, Aa it01 the corresponding statistical errors due to numerical 
integration. This allows to quantify the mass shift due to 0{a) corrections. The shift due to 
higher-order corrections is derived according to the same procedure, by generating a sample 
of pseudo-data for the Mt distribution at 0(a) and fitting them in terms of the Mr spectrum 
obtained including higher-order corrections for 10 hypothesized W mass values. In this case, 
we use 1 MeV spacing between masses. Figure E] shows the A% 2 = x 2 ~ Xmin distributions 
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FIG. 2: The Ax 2 = X 2 — Xmin distributions from a fit to the Mt distribution, including 0(a) QED 
corrections (left) and higher-order QED corrections (right), as a function of the W mass shift, at 
yfs = 2 TeV. The results for the W — > ev and W — > /U^ channels are shown. 

as a function of AMw = — , for the fit with 0(ai) corrections (left) and the fit 
with higher-order corrections (right). The mass shift observed for 0(a) corrections amounts 
to about 20 MeV for the W — > ev decay (dashed line) and to 110 MeV for the W — > fiv 
decay (solid line), as a consequence of the different identification requirements. These shifts 
are in reasonable agreement with the results of the CDF and D0 collaborations, even in 
the absence of a complete detector simulation. The mass shift due to higher-order effects 
is about 10 MeV for the W — > \xv channel (solid line) and a few MeV (dashed line) for the 
W — > ev channel. We performed the same analysis for the LHC collider (using the cuts and 
pseudo-detector simulation of the Tevatron collider) and found that the same conclusions 
do apply to the LHC. 

In conclusion, we have evaluated the impact of higher-order final-state QED corrections 
on the determination of the W mass at hadron colliders, in view of future improved measure- 
ments with an accuracy of 15-30 MeV. In the presence of realistic selection criteria, we have 
found that the shift due to these corrections is about 10 MeV in the W \iv channel and 
practically negligible in the W — > ev channel. The calculation, if included in future experi- 
mental analyses, would reduce the uncertainty in the precision measurement of the W mass 
at hadron colliders. To this end, the Monte Carlo program HORACE is available for data 
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analysis. A more realistic analysis would require a full detector simulation, which is beyond 
the scope of the present paper. The study of the neutral-current process pp — > 7, Z — ► l + l~ 
is left to a future work. 
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